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Abstract 

In  this  paper  we  study  the  response  to  external  disturbances  of  a  double  integrator  with  saturating  feedback.  For  a  class  of  disturbances 
that  have  bounded  integrals  over  all  intervals,  we  show  that  a  linear  static  feedback  law  can  always  be  designed  to  ensure  boundedness 
of  the  states.  Moreover,  boundedness  can  be  preserved  if  the  disturbance  is  biased  by  a  small  DC  signal.  In  the  special  case  that  the 
disturbance  is  made  up  of  a  finite  number  of  sinusoids  and  a  small  bias,  any  linear  static  feedback  ensures  boundedness  of  the  states. 
These  results  are  an  extension  of  previous  results  by  Wen,  Roy,  and  Saberi  (2008). 
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1  Introduction 

The  double  integrator  system  is  commonly  seen  in  control 
applications  including  low-friction,  free  rigid-body  motion, 
such  as  single-axis  spacecraft  rotation  and  rotary  crane  mo¬ 
tion  (see  Rao  and  Bernstein,  2001,  and  references  therein). 
Of  particular  interest  is  the  control  of  double  integrators 
subject  to  input  saturation.  A  classical  result  is  that  a  dou¬ 
ble  integrator  with  a  saturating  linear  static  feedback  pro¬ 
vides  global  asymptotic  stability  of  the  origin.  This  result 
has  been  extended  to  mixed-type  systems  by  Tyan  and  Bern¬ 
stein  (1999)  and  Yang,  Stoorvogel,  and  Saberi  (2010).  Many 
other  control  methods  have  also  been  proposed.  A  brief  sum¬ 
mary  and  comparison  of  various  methods  is  given  by  Rao 
and  Bernstein  (2001). 

Compared  with  the  relatively  mature  study  of  internal  sta¬ 
bilization,  the  dynamic  response  of  a  double  integrator  with 
saturating  feedback  to  external  disturbances  is  still  not  fully 
understood.  In  this  paper,  we  study  the  disturbance  response 
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of  a  double  integrator  controlled  by  a  saturating  linear  static 
state  feedback,  as  given  below: 


X2  =  G  {—k\X\  —  k,2Xl)  +  d . 

where  a  represents  the  standard  saturation  function  o(s)  = 
s/gn(s)min{l,  |s|}.  The  goal  is  to  identify  a  class  of  distur¬ 
bances  for  which  the  states  of  the  above  controlled  system 
remain  bounded.  Chitour  (2001)  and  Stoorvogel,  Shi,  and 
Saberi  (2004)  have  previously  studied  this  problem  in  the 
context  of  f£p  stability.  They  showed  that  for  any  k  \  >  0  and 
ko  >  0,  (1)  is  2z?p  stable  for  p  G  [1,2],  and,  moreover,  the 
trajectories  remain  bounded  for  any  d  G  f£p  with  p  G  [1,2]. 
However,  for  any  p  >  2  there  exists  d  G  2z?p  that  cause  the 
states  to  grow  unbounded  from  certain  initial  conditions.  Shi 
and  Saberi  (2002)  studied  the  notion  of  input-to-state  sta¬ 
bility  (ISS)  (see  Sontag  and  Wang,  1995)  for  system  (1), 
and  proved  that  no  choice  of  k\,k2  can  achieve  ISS.  Specifi¬ 
cally,  there  exist  bounded  disturbances  with  arbitrarily  small 
Jzfoo  norm  that  cause  the  states  to  grow  unbounded  from  cer¬ 
tain  initial  conditions.  Even  more  dramatically,  unbounded 
growth  can  be  achieved  by  vanishing  disturbances  with  ar¬ 
bitrarily  small  Jzfoo  norm. 

Wen  et  al.  (2008)  extended  the  negative  result  from  Shi  and 
Saberi  (2002)  to  classes  of  small  JzfL  signals  with  further 
restrictions.  However,  Wen  et  al.  (2008)  also  showed  that 
boundedness  of  the  states  of  (1)  is  preserved  for  a  particular 
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class  of  small  disturbances  that  have  bounded  integrals  over 
all  intervals.  To  be  precise,  let  a  family  of  parameterized  sets 
be  defined  by 


Proof.  Define  y i  =  k\X\  -\-k2X2,  yi  =  &2T2  and  t  =  ^t.  Then 
the  closed-loop  system  in  the  new  coordinates  and  with  t  as 
the  time  variable  becomes 


CIm  = 


j*/€  JSf«(l) 


Vfi,f2  >  0, 


where  jSfL>(l)  denotes  the  set  of  signals  of  magnitude 
less  than  1 .  For  a  given  linear  state  feedback,  it  was  proven 
that  for  any  M  >  0  there  exists  a  q*  >  0  such  that  the  states  of 
(1)  remain  bounded  for  any  disturbance  d{t)  =  qs(t ),  where 
s(t)  £  Q.m  and  q  <  q*.  In  other  words,  signals  with  bounded 
integrals  can  be  tolerated  if  they  are  scaled  down  by  a  suf¬ 
ficient  amount.  It  was  furthermore  shown  that  a  small  bias 
in  d(t)  can  be  tolerated  while  still  achieving  boundedness 
of  the  states.  The  class  of  disturbances  considered  by  Wen 
et  al.  (2008)  covers  a  broad  class  of  signals,  such  as  peri¬ 
odic,  quasi-periodic,  and  .2)  signals. 


-Zl(t)=y2(t)-X[a(yi(t))-d(t)}, 

d  (3) 

where  A  =  k2/k\  >0.  We  shall  construct  an  integral- 
bounded  disturbance  d  that  causes  the  states  to  grow 
unbounded  from  a  particular  initial  condition. 


Step  1:  Suppose  the  trajectory  of  (3)  starts  from 
A  =  (l, NX),  for  some  large  integer  N,  at  time  tx  =  0. 
We  will  construct  a  d  to  drive  the  states  from  point  A  to 
B  =  (1,  — (AA  +  2)A)  at  time  ?g.  Choose 


This  paper  is  an  extension  of  the  work  in  Wen  et  al.  (2008). 
Our  focus  is  also  on  disturbances  with  bounded  integrals; 
however,  we  strengthen  the  results  from  Wen  et  al.  (2008)  by 
dispensing  with  the  magnitude  restriction  on  d  and  removing 
the  attenuation  factor  q.  Specifically,  we  consider  the  family 
of  parameterized  sets 


klM  — 


|c/  £  J£oo  |  Vfj  ,?2  >  0, 


<  M  >  . 


(2) 


d(t)  =  2nsm(m). 


d\(t)  =  [  d(r)dx  =  2(1  —  cosOrf)), 
Jo 

d2(t)  =  [  di(r)dT  =  2t—  |  sin(7rf). 
Jo 


If  a  signal  d  belongs  to  CLm  for  some  M  >  0,  we  refer  to 
it  as  an  integral-bounded  signal,  and  we  refer  to  M  as  an 
integral  bound  on  d. 

We  first  prove  a  new  negative  result,  namely,  that  for  a 
given  linear  static  feedback  law,  there  always  exist  integral- 
bounded  signals  that  cause  the  trajectories  of  the  system  to 
grow  unbounded  from  some  initial  conditions.  Our  next  re¬ 
sult,  however,  shows  that  if  an  integral  bound  M  is  known 
a  priori,  then  k\  and  kj  can  always  be  designed  to  ensure 
boundedness  of  the  states,  regardless  of  initial  conditions. 
Moreover,  boundedness  can  be  ensured  also  if  the  integral- 
bounded  disturbance  is  biased  by  a  DC  signal  of  magnitude 
less  than  1 .  Finally,  we  prove  an  even  stronger  result  for  dis¬ 
turbances  consisting  of  a  finite  number  of  sinusoids  plus  a 
DC  offset  of  magnitude  less  than  1.  In  this  case,  any  inter¬ 
nally  stabilizing  linear  static  feedback  law  ensures  bound¬ 
edness  of  the  states. 

2  Main  result 

The  first  theorem  shows  that  not  every  internally  stabilizing 
static  linear  law  can  maintain  boundedness  of  the  trajectories 
in  the  face  of  integral-bounded  disturbances. 

Theorem  1  Consider  the  system  (1)  with  A 1  >  ( )  and  A  '  >0. 
There  exists  an  integral-bounded  signal  d  and  an  initial 
condition  such  that  X\  and  X2  grow  unbounded. 


Since  d,vA  (0)  =  3^2(0)  —  A,  we  see  that  for  large  N,  the  tra¬ 
jectory  will  initially  move  to  the  right.  If  y\ (t)  >  1,  we  have 

22(1)  =  3^2(0)  —  Xt  +  Xd\{t)  =NX  —  Xt  +  2A(1  —cos  (nt)) 

and 

yi (t)  =  yi (0)  +  /  y2{r)dt-Xt  +  Xdi(t) 

Jo 

=  yi (0)  +NXt-  \t2+Xd2{t) -Xt  +  Xdi(t) 

=  yi  (0)  +NXt—  \r  +  2Xt-  sin(zrr)  —  Ar 

+  2A(1  —cos  (nt)) 

=  —  jt2  +  (NX  +  X)t  +  2A(1  —cos  (nt)) 

-  ^  sin(7K)  +  1. 

Given  a  sufficiently  large  N,  y\ (?)  only  has  one  intersection 
with  vi  =  1  for  t  >  0.  This  is  shown  in  the  Appendix.  For 
t  =  2N  +  2,  we  have 

yi  =-X(1Nf-)2  +(AA+A)  x  (2A  +  2)  +  1  =  1 

and 

r2N+2 

yi  =  NX~X{2N  +  2)+X  /  d(t)dt  =  -(N  +  2)X. 

Jo 
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This  shows  that  the  trajectory  will  cross  y\  =  1  at  B  = 
(1 ,  —  (N  +  2)A)  at  time  ?g  =  2N  +  2.  We  have 


<  4  with  tA<h<h<  tB- 


Step  2: 


From  fg,  set  d  (?) 

dyi 

d? 


=  Ct(yi(f)).  Then 

(t)=y  2, 


dy2 

d t 


{i)  =  0. 


The  trajectory  will  move  directly  toward  the  left  to  C  = 
(—  1,—  (iV  +  2)A)  at  time  tc  =  ?g  +  <NA2)X  •  Clearly,  for  fV 
sufficiently  large, 


<  1  with  ?g  <  ?i  <  ?2  <  ?c- 


y2 


Step  3:  From  ?c,  choose 


Fig.  1.  State  trajectories  of  double  integrator  with  A  =  1 


= 2tt sin f 7tF) .  for  any  0  =  tA  <  ?i  <  ?2-  In  the  original  time  variable  t,  we 

calculate  the  integral  bound  as 

Following  the  same  argument  as  in  Step  1  with  N  replaced 
by  N  =  N  +  2,  we  find  that  the  trajectory  will  re-cross  y\  = 

—  1  at  D  =  (—  1,  (iV  +  2)A)  =  (—  1,  (IV  +  4)A)  at  time  to  = 
tc  +  2N  +  2  =tc  +  2N  +  6.  Similarly 


<  4  with  ?c<?i  <  ?2  <  to  ■ 


Hence,  we  have  shown  that  there  exists  an  integral-bounded 
signal  with  integral  bound  M  =  M  =10^  that  causes  the 
states  to  grow  unbounded  for  a  particular  initial  condition.  ■ 


Step  4:  From  to,  choose 

d(t)=(j{yi{t)). 

The  trajectory  will  move  directly  to  the  right  and  cross  y i  =  1 
at  E  =  (1,  (IV +  4) A)  at  time  ?£=?£>+  jfyy-  We  have 


Next,  we  show  that  if  an  integral  bound  M  is  given  a  priori , 
then  we  can  always  design  a  static  stabilizing  linear  feedback 
to  ensure  boundedness  of  the  trajectories. 

Theorem  2  Let  M  be  given.  Ifk\  and  k2  satisfy  ^  >  16 M, 
then  for  any  d  €  LIm  and  any  initial  condition,  we  have 
x\.x2  € 


<  1  with  to  <  ?i  <  ?2  <  ?£• 


Proof.  The  proof  of  Theorem  2  is  a  consequence  of  Lemmas 
3  and  4  which  are  stated  and  proved  below. 


.  Lemma  3  Consider  the  system 
The  system  trajectory  resulting  from  Steps  1  through  4  is  vi¬ 
sualized  in  Fig.  1.  By  repeating  these  steps,  the  state  grows 

unbounded,  and  we  can  check  that  the  constructed  distur-  X|  —  x2+y,  ^ 

bance  signal  satisfies  *2  =  cr(— Ai^ri  —k2x2), 


where  |y(f)|  <  2 M  for  all  t  >  0  and  ^  >  16M.  In  that  case, 
we  have  x\  ,x2  £  Jf^for  any  initial  condition. 
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Proof  of  Lemma  3.  Define  a  positive  definite  function  V  as 


Lemma  4  Consider  the  system 


This  function  was  first  introduced  in  Chitour  (2001).  Differ¬ 
entiating  V  along  the  trajectories  yields 


V  =  (k\x2  +  k\y)c  (k\X\)  —  2k\x2o(k\x\  +  k2x2) 

+  [k\X2  +  k\y  -  k2o(k\x\  +  k2x  2)]  o(k\x\  +  k2x2) 

=  k\X2  \o[kix\)  —  a{k\xi+k2X2)\—k202{k\x\  +  k2x 2) 

+  £iy[(7(£ixi  +  &2x2)  +  c(fcixi)] 

<  &ix2  [ct(Atixi)  -  <7(&ixi  +  &2x2)] —  k2o2{k\x\  +  k2x 2) 

+  2&i|y|. 


If  |£ixi  +fc2x2|  >  5,  then 

— ^2Ct2(^ixi  +£2X2)  +  2£i|y|  <  —  16Mk\  x  \  +  4k\M  <  0. 
Hence 


y  =  o{v)-c(y+k2y)  +  d,  y(0)  =  0,  (5) 

where  k2  >  0,  d  G  LIm  and  v  is  continuous.  We  have  |y(f)|  < 
2 M  for  all  t  >  0. 

Proof  of  Lemma  4.  Define 

y  =  d,  y(0)  =  0. 

Since  d  G  the  solution  satisfies  |y(f)|  <  Af  for  all  1  >  0. 
Define  y  =  y—y.  We  have 

^  =  a(v)-c7(v+fc2(y+y)),  y(0)  =  0. 

Define  a  positive  definite  function  V  =y2.  Taking  the  deriva¬ 
tive  of  V  with  respect  to  t,  we  get 

V  =  9  [o  ( v)  -  a  ( v  +  k2  (y + y) )] . 


V  <  k\X2  [ct(^ixi)  —  o{k\x\  +  k2x2)\  <  0. 

If  I&1X1+&2X2I  <2  ,  then  by  using  Lemma  B.2  in  Shi,  Saberi, 
and  Stoorvogel  (2003),  we  get 

&1X2  [<j(kiXl)~  o{k\X\  +£2X2)]  <  -yX2Cj(&2X2). 

If  we  also  have  that  |x2|  >  max{8M,  ^},  then 

k\x2  [(7(feixi)  —  a{k\x\  +  fc2x2)] 

<  -yX2Cr(^2X2)  <  —4k\M, 

which  yields  V  <  0.  We  therefore  conclude  that  V  <  0  out¬ 
side  the  region  defined  by  \k\X\  +fc2x2|  <  \  and  |x2|  < 
max{8M,  ^}.  Hence,  V  remains  bounded,  which  implies 
that  xi,x2  €  ■ 

Now  consider  the  double  integrator  system  (1).  We  construct 
a  fictitious  state 

y  =  o{-k\xi  - fc2x2)  -  cr(  k\x\  - &2x2  +k2y)  +d, 

y(  o)  =  o. 

By  defining  z  =  x2  —  y,  we  obtain  the  augmented  system 

xi  =z+y, 

z  =o{—k\X\  -k2z), 

y  =c(—k\x\  -k2x2)-  a{-k\x\  -  k2x2  +  k2y)  +  d , 

with  y(0)  =  0,  z(0)  =  X2(0).  From  Lemma  3,  we  know  that 
given  >  16M,  x\  and  z  remain  bounded  provided  |y|  <2 M. 
The  latter  statement  is  proven  by  the  following  lemma. 


IfV>M2,then  \y(t)\  >M>  |y(f)|,  which  implies  that  k2  (y  + 
y)  has  the  same  sign  as  v.  It  then  follows  that  V  <  0.  Since 
V (0)  =  0,  we  can  conclude  that  V  <  M 2  and  thus  ]y(f)|  <  M 
for  all  t  >  0,  and  it  follows  that  |y(f)|  <  |y(f)|  +  |y(f)|  <  2 M 
for  all  t  >0.  m 

From  Lemmas  3  and  4,  we  know  that  x\  and  2  are  bounded. 
Since  y  is  bounded  as  shown  in  Lemma  4,  we  conclude  that 

xi,x2  G  _S?L>.  ■ 

An  immediate  consequence  of  Theorems  1  and  2  is  that  if  k\ 
and  k2  are  arbitrary  positive  real  numbers,  then  boundedness 
is  guaranteed  if  the  integral  bound  M  is  sufficiently  small. 
This  is  formally  stated  in  the  following  corollary. 

Corollary  5  For  any  given  k\  >  0  and  k2  >  0,  we  have 
xi,x2  G  J2L  if  d  G  CIm  with  M  <  j^-. 

In  the  next  theorem  we  consider  integral-bounded  distur¬ 
bances  that  are  biased  by  a  DC  signal.  We  show  that,  if  the 
magnitude  of  the  bias  is  less  than  1  by  a  known  margin,  and 
an  integral  bound  M  is  known  a  priori ,  then  k\ ,  k2  can  be 
chosen  to  ensure  boundedness  of  x\  ,x2. 

Theorem  6  Let  M  >  0  and  8  G  (0, 1]  be  given,  and  suppose 
that  d  =  d\  +  d2  where  d\  is  a  constant  with  |t/i|  <  1  —  8 
and  d2  G  LIm-  If  k\ ,  k2  satisfy  k2  >  max{  ,  48gj>M  },  then 
xi,x2  G  J2L. 

Proof.  The  closed-loop  system  is  given  by 

xi  =x2, 

x2  =  a{—k\xi  -  k2x2)  +d\  +d2. 
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We  construct  a  fictitious  state 

y  =  o(—k\xi  -  k2X 2)  -  o(-k\x\  - £2*2  +  kzy)  + tfe, 

y(  o)  =  o. 

Lemma  3  shows  that  |y(f )  |  <  2M  for  all  t  >  0.  Similar  to  the 
proof  of  Theorem  2,  we  define  z  =  x 2  —  y  and  convert  the 
closed-loop  system  to  the  form 

x\  =  z+y, 

z  =  c{—kixi  —  k2z)  +d  1 , 
with  z(0)  =  X2  (0)  and  >’  £  Jz?oo(2 M). 

We  also  introduce  another  fictitious  state 

vv  =  cr(— Ari-Xi  —k2z)  —  cr(— £1X1  —  &2Z  +  &2W-  ^1) 


(3)  s  [o-dl  (v  +  s)  -  &dl  (v)]  >  0 

Moreover,  it  is  shown  in  Lemma  9  in  the  Appendix  that  if 
|v|  <  |,  then 

s  [odi  (v  +  s)~  &di  (v)]  >  sgS/2(s) , 

where  g§/2  (s)  is  the  standard  saturation  function  with 
saturation  level  8/2,  which  is  defined  by  C7s/2(S)  — 
sign(s) min{(5/2,  |s|}. 

With  this  generalized  saturation  function,  the  closed-loop 
system  can  be  rewritten  as 

4 1  =  fe+w  +  y, 

4 2  =  Gdi{-ki£i  —k2%  2). 


with  w(0)  =  0.  Following  the  same  argument  as  in  the  proof 
of  Lemma  4,  we  can  show  that  |w(/)|  <  <  M  for  all 

f  >  0.  Define  =x\,  ^2  =z  —  n’  =  x2  — y  —  w.  Then  (1)  can 
be  transformed  into 


Define  a  positive  definite  function 


Gdx 


41  =42  +  H’  +  y, 

42  =  <?(— &i4i  —  fe42 — ^1)  +^t, 

where  <jji(0)  =  xi(0),  42(0)  =  x2(0)  and  |w(f) +y(f)|  <M  + 
2 M  =  3 M  for  all  t  >  0.  Since  w(f)  and  y(t)  are  bounded,  we 
know  that  x\  and  x2  are  bounded  if  <*]  and  5?  are  bounded. 

Define  adl  (s)  =  c(s  —  d\)+d\  with  \d\  \  <  1  —  8.  Then 

{1  d \ .  s  ^  1  -(-  d\ , 
s,  -1  +di  <  s  <  1  +di,  (6) 

—  1  -{ ~ d\,  s  ^  —  l+d\. 

This  function  can  be  viewed  as  a  generalized  saturation  func¬ 
tion,  which  is  visualized  in  Fig.  2.  It  is  easy  to  verify  that 


Differentiating  V  along  the  trajectory  yields 

V  =  [k\^2+k\w  +  k\y)acjl  (fci4t)  -  2ki&cdl  (*t4t  +^42)+ 
\k\t,2+k\w  +  k\y  —  k2ddl (&i4i  +  k2£i)\  Gdl  (^i4i  +^42) 
=  k\t,2  \?di  (^i4i)  —  oi/i  (*t4t  +^242)]  -k2ajl  {k\E,\  +  £24 2) 

+  £i(w+y)  [5-dl (^t^O  +  Si/j (&i4i  +^242)] 

<  ^1^2  [5-rfl  (*i4i)  —  a-dl(A:i4i  +^2^2)]  —k2Gdi(ki£  1  +^42) 

+  12&iM. 


If  |^i %\  +  £24 2I  >  f,  then 


-*28*  (*i|i  +  ^2)  +  12*iM  <  -  ^  f  +  \2k\M  =  0, 

and  hence  V  <  0.  If  |&i4i  -F  ^2 42 1  <  §  and  |42|  > 
max{  2f,2 , }.  then  by  using  Lemma  9  we  have 

£i42  [8rfl(fci4i)-arfl (&i4i  +  £24 2)]  < 

-k&Gspfab)  <-ki^<  —Y2k\M, 

and  hence  V  <  0.  We  therefore  find  that  V  <  0  outside  the  re¬ 
gion  defined  by  \kl%l+k2&\  <  \  and  |42|  <  max{^ ,  ^}. 
It  follows  that  V  remains  bounded,  which  implies  that  4i 
and  Q  remain  bounded.  ■ 

Our  final  result  concerns  a  special  case  where  the  disturbance 
consists  of  a  finite  number  of  sinusoids  together  with  a  DC 
bias  of  magnitude  less  than  1.  In  this  case,  any  internally 
stabilizing  linear  static  feedback  controller  guarantees  that 
the  states  of  the  system  (1)  remain  bounded. 


Gdx  satisfies  the  following  properties: 

(1)  |a*Cs)|<2 

(2)  sdd j  (s)  >  0  and  sGdl  (s)  =  0  iff  s  =  0 
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Theorem  1  Consider  the  system  (1)  with  k\  0  and  It2  >0. 
Suppose  that  d  =  d\  +  d2,  where  d\  is  a  constant  satisfying 
|c/l|  <  1  and  d2  is  generated  by  an  exogenous  system 


Appendix 

Intersection  problem  in  the  proof  of  Theorem  1 


vv=Aw,  w(0)=wo, 
d  =  Cw, 


We  shall  show  that 


where  A  is  non-singular  and  satisfies  A  +A'  =  0.  We  have 
Xi,X2  €  Afoofor  any  initial  condition. 

Proof.  We  can  rewrite  the  closed-loop  system  in  a  compact 
form: 


vv 

A  0  0 

w 

0 

x\ 

= 

0  0  1 

xi 

+ 

0 

X2_ 

COO 

*2_ 

1 

[a(-k\x\  —k2X2)  +d\] . 


yi(f)  =  —  ^t2 +  (NX+X)t+2X(l  —  cos(nt))  —  ^  sin(7rf)  +  l 

has  only  one  intersection  with  y  i  =  1  for  t  >  0  and  suffi¬ 
ciently  large  N.  Let 


t\  =  min{f  >  0 : yi (f)  =  1},  Fq  =  min  ^  t  >  0  :  (t )  =  0  ^  . 


Given  vi  (0)  =  1  and  ^  >  0,  we  must  have  t\  >Iq.  Note  that 


Consider  the  state  transformation 


w 

1 

o 

o 

1 _ 

vv 

xi 

= 

~CA~2  1  0 

X\ 

X2_ 

—CA~1  0  1 

X2_ 

This  transformation  results  in  the  system 


vv 

A  0  0 

w 

0 

x\ 

= 

0  0  1 

xi 

+ 

0 

_X2_ 

0  0  0 

x2_ 

1 

x  [a  {-{kxCA~2  +  k2CA~l)w  -  km  —  k2x2)  +di]  . 

Define  v=  —  (k\CA~2 -\-k2CA~x)w +d\.  Then 

o(—(k\CA~2  +k2CA~l)w  —  kixi  —k2x2)  +  d\ 

=  c(v  —  k\x\  —  k2x2  —  d\)  +d\ 

=  &dl  {—k\x\  -  k2x2  +  v), 

where  adl  is  the  generalized  saturation  function  defined  in 
the  proof  of  Theorem  6.  The  dynamics  of  xi  and  x2  can  now 
be  written  as 


=  —  Xt+NX  +  X+2Xnsin(nf)  —  2Xcos(jtt). 

Hence  to  >  N  +  1  —  2n  —  2  >  j  for  sufficiently  large  N. 
However, 

y\ (t)  >  -  § 12  +  (NX  +  X)t- 4A  - 2A  +  1 
Hence  we  have 

—  +  (NX  +  A )F i  —  4A  —  2A  <0 

or  equivalently 

^'r-(A+l)Fi+6>0 

This  implies 

h<n,  or  t\  >  r2 

where  r\  and  r2  are  two  roots  of  ^t2  —  (N  +  l)f  +  6. 
n,2=N+lT^J(N+l)2-l2 

Note  that 


xi  =x2, 

X2  =  odl  (-k\X\  -  k2X2  +  v). 

Clearly  v  €  JKo.  It  was  shown  by  Chitour  (2001)  that  the 
(xi  ,x2)  dynamics  is  0  stable  from  v  to  Jci  and  x2  for  any 
k\  >  0  and  k2  >  0.  ■ 

Remark  8  For  ease  of  presentation,  we  use  a  standard  sat¬ 
uration  function  with  saturation  level  1,  but  all  the  results 
obtained  in  this  paper  can  easily  be  extended  to  the  case 
where  a  saturation  function  with  arbitrary  saturation  level 
A  is  used. 


n  =n+i- yV+i)2  — 12 

12  N 

N+l  +  s/(N+l)2~l2  2 

for  sufficiently  large  N.  Since  we  already  know  h  >  i0  >  f , 
we  must  have  fi  >  r2.  But  then 

r2=N+l  +  ^J(N+l)2-l2>N+l  +  \/(N+l)2/4 

=  l(N+  1) 


6 


for  large  N.  We  find  that  t\  >  |  (N  + 1) .  But  for  t  >  |  (N  + 1 ), 
we  have 

^L(f)  <  —  jX(N+  1)  +NX+X  +2A?rsin(?rf)  —  2Xcos(nt) 
<  —\X(N+  \)+2Xnsm(nt)  — 2Acos (nt)  <  0 

for  sufficiently  large  N.  This  shows  that  y\  (t)  <  1  for  all 
t>t\,  and  hence,  the  only  intersection  with  vi  =  1  is  at  t  t\ . 

Property  of  6d]  (s) 

Lemma  9  The  generalized  saturation  function  dd{  defined 
in  (6)  with  \d\  \  <  1  —  5  satisfies 

s  [ad j  (s  +  v)  -  6dl  (v)]  >  so5/ 2  (s) 

for  |v|  <  j  where  <75/2(5)  denotes  the  standard  saturation 
function  with  saturation  level  8/2  defined  as  <75/2(5)  = 
sign(s)  min{5/2,  |5|}. 

Proof.  If  |s|  <  we  have  |v+5|  <  8  <  1  —  \d\  |.  By  definition 
(6),  we  have  ad{  (5  +  v)  =  5  +  v.  Hence 

adl  (5  +  v)  -  ad j  (v)  =  5  +  v-  v  =  s. 

If  |s|  >  f ,  it  can  be  seen  from  Fig.  2  that 

I drfj  (5  +  v)  —  adl  (v) |  >  \sign(s)j  +v  — v|  =  f. 

Hence  5  [<7dl  (5  +  v)  -  adl  (v)]  >  sas/2(s).  m 
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